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Abstract 

The Luther-Emery liquid is a state of matter that is predicted to occur in one-dimensional 
systems of interacting fermions and is characterized by a gapless charge spectrum and a gapped 
spin spectrum. In this Letter we discuss a realization of the Luther-Emery phase in a trapped 
cold-atom gas. We study by means of the density- matrix renormalization-group technique a two- 
component atomic Fermi gas with attractive interactions subject to parabolic trapping inside an 
optical lattice. We demonstrate how this system exhibits compound phases characterized by the 
coexistence of spin pairing and atomic-density waves. A smooth crossover occurs with increasing 
magnitude of the atom-atom attraction to a state in which tightly bound spin-singlet dimers 
occupy the center of the trap. The existence of atomic-density waves could be detected in the 
elastic contribution to the light-scattering diffraction pattern. 
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Introduction — The ability to confine ultracold Bose and Fermi gases inside artificial 
crystals generated by standing- wave laser light fields, i.e. optical lattices (OLs), offers the 
possibility to create ideally clean and highly tunable strongly interacting quantum many- 
body systems Q]. The experimental observation by Greiner et al. Q] of a superfluid-Mott 
insulator transition in a 3D OL loaded with 87 Rb atoms demonstrated the possibility of 
realizing strongly correlated bosonic systems. A number of experimental results have also 
been obtained for fermionic systems in OLs jjj. Cold atoms have been successfully trapped 
in low-dimensional geometries. A 87 Rb gas has been used to realize experimentally a Tonks- 
Girardeau gas U The preparation of two-component Fermi gases in a quasi-l-D geometry [^] 
provides a unique possibility to experimentally study phenomena predicted a long time ago 
for electrons in ID. Spin-charge separation in Luttinger liquids is a paradigmatic example p. 

In the presence of attractive interactions fermions are predicted to form a peculiar ID 
liquid phase characterized by a massive spin sector, i.e. a Luther- Emery liquid [7|. The 
gap in the spin sector induce, an exponential decay of sp.n conations, while sing.et su- 
perconducting and charge- density wave correlations have a power-law decay So far no 
observation of the Luther-Emery phase has been reported in solid-state electronic systems. 
In Ref. [8] it has been shown that an integrable model of two-component interacting Fermi 
gases in a quasi-l-D geometry exhibits a smooth crossover between a Luther-Emery liquid 
and a Luttinger liquid of tightly-bound bosonic dimers. Although long-range phase coher- 
ence is absent in a ID Luther-Emery liquid, antiparallel-spin fermions do pair as they do in 



a conventional superconductor. Seidel and Lee jgj have shown that the ground-state enertiy 
of a system with spin gap and gapless charge degrees of freedom has an exact period of 
he/ (2e) (corresponding to half a flux quantum) as a function of an applied Aharonov-Bohm 
flux. 

In this Letter we propose that ID Fermi gases in OLs be used to study the Luther- 
Emery phase and the existence of antiparallel-spin pairing in a Luther-Emery liquid. A 
realistic OL is superposed on a slowly varying harmonic potential that makes the lattice 
sites inequivalent. In order to assess the very existence and the nature of the Luther- Emery 
phase in realistic cold-atom systems it is thus of fundamental importance to analyze the 
interplay between attractive interactions and the confining potential. This is not a merely 
quantitative issue: it is well known that in other cases the confinement modifies qualitatively 
the properties of the gas. An example is the coexistence of superfluid (metallic) and Mott- 
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insulating regions in bosonic (fermionic) OLs 

We first show that, in the presence of harmonic confinement, a ID Fermi gas with attrac- 
tive interactions inside an OL manifests unambiguous real-space spin pairing, which in turn 
determines the emergence of Atomic-Density Waves (ADWs) in the ground-state density 
profile. We then propose an experiment that can lead to the observation of these ADWs. 

The inhomogeneous Fermi-Hubbard model — We consider a two-component Fermi gas 
with N atoms confined by a harmonic potential of strength V2 in a ID lattice with unit 
lattice constant and L lattice sites i G [1,L\. The system is described by an inhomogeneous 
Fermi-Hubbard Hamiltonian, 

n = -t ^(4c i+lff + h.c.) + uJ2 + ^ - L / 2 ) 2 ^ • (!) 

i,a i i 

Here t is the hopping parameter, a =j, J, is a pseudospin-1/2 degree of freedom (hyperfine- 
state label), U is the on-site Hubbard interaction parameter, hi = ^2 a hia — X] ff C^ anc ^ 
hi) = N. L is always chosen so that the trap makes the ground-state site occupation 
n i — (^Gsl^il^cs) go to zero smoothly near the edges of the lattice. The Hamiltonian 7i with 
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repulsive interactions (U > 0) has been extensively discussed in the literature 
In this Letter we focus on attractive interactions, i.e. U < 0, and on the unpolarized case 
(AT T = N^. 

In the absence of harmonic confinement (V2 = 0) TC has been exactly solved by Lieb and 
Wu 14] by means of the Bethe Ansatz and shown to belong to the Luther-Emery liquid 
universality class. For V2 7^ we calculate the ground-state properties of H in Eq. ((TJ) by 
resorting to the density-matrix renormalization-group (DMRG) method [15], which provides 
a practically exact solution for any value of U/t. Motivated by the recent interest in the 
dev elop ment of density-functional schemes for strongly correlated ID systems jl^ . IlI llfil 
Q Q- Q we also use, in parallel to DMRG, a lattice density-functional scheme based on 
the Lieb-Wu solution for V2 = 0. These calculations employ the Bethe-^nsafe Local-Density 
Approximation (BALDA) in its fully numerical formulation Q] (for odd iV the ground 
state has a finite polarization and thus a generalization to spin-BALDA is necessary). The 
BALDA method allows us to efficiently treat complex systems with a very large number 
of sites. The combined use of DMRG and BALDA allows a detailed understanding of the 
problem and paves the way for further extensions. 

The pairing gap — We first discuss the tendency to pairing by analyzing the pair binding 



TABLE I: Ground-state and pair-binding energies for N = 30, L = 100, and Vi/t = 4 x 10 -3 . 
The agreement between DMRG and BALDA for E G $ is quite satisfactory even for U/t = —20, 
where the deviation is about 2.2%. However, BALDA tends to overestimate E-p even at moderate 
values of U/t. The "x" sign indicates that the spin-BALDA code for 31 atoms does not converge 
for U/t = -20. 



U/t E»f OA /(tL) E™ RG /(tL) E* ALDA /t E° MRG /t 


-0.5 -0.35824 


-0.35832 


-0.0283 


-0.0213 


-1 -0.39336 


-0.39340 


-0.0614 


-0.0577 


-2 -0.47672 


-0.47631 


-0.3265 


-0.2442 


-4 -0.70693 


-0.69010 


-5.1008 


-1.3278 


-20 -3.05320 


-2.98536 


x 


-16.4217 



binding energy of finite 
In Table|Uwe report 



energy defined as 

E P = E GS (N + 2) + E GS (N) - 2 E GS (N + 1) (2) 

where E GS = (^gsI^I^gs) is the ground-state energy. The pair 
Hubbard chains with attractive interactions has been studied in Ref . _ 
results for A^ = 30 fermions in a lattice with L = 100 sites, inside a trap with V^/t = 4x 10~ 3 . 
E-p is negative, signaling a tendency to pairing. 

Atomic- density waves — The pairing is associated with the presence of ADWs, which are 
stabilized by the harmonic potential. In Fig. ^ we report our numerical results for the site 
occupation of a gas with A^ = 30 atoms. The consequences of Luther- Emery pairing in the 
presence of confinement are dramatic. For U < the site occupation exhibits a density wave 
(with Af/ 2 peaks in a weak trap), reflecting the tendency of atoms with different pseudospins 
to form stable spin-singlet dimers that are delocalized over the lattice. 

For small Vijt (see the top panel of Fig. ^) in the bulk of the trap (80 < i < 100) 
can be fitted to an ADW of the form rij = n + A ADW cos (lc ADW i + if). For example, for 
V 2 /t = 10~ 5 we find k ADW = 0.73 for U/t = -1 and k ADW = 0.84 for U/t = -3. In 
such a weak confinement the oscillations of the site occupation extend to regions far away 
from the center of the trap, where they are characterized by smaller edge wavenumbers. 

n 

For V2 = bosonization predicts |6| an incipient ADW with wavenumber 2kp, kp being 



the Fermi wavenumber, quenched by strong quantum fluctuations. In the present case we 
find A; A dw — 2fcp ff , where the effective Fermi wavenumber fc F ff = 7rn/2 is determined by the 
average density in the bulk of the trap (note that k ADW = n when the average density in the 
bulk reaches half filling). 

Finite-size effects become important on increasing V2 (see the bottom panel of Fig. 
and a simple fitting formula such as the one used above does not work even at the center 
of the trap. Eventually when Vijt ~ 10 _1 a region of doubly-occupied sites develops at the 
center of the trap: spin-singlet dimers, which in a weak trap are delocalized, are squeezed 
close together to produce an extended region of ~ N/2 doubly-occupied consecutive sites. 

In Fig. |2]we show how the ADWs evolve with increasing \U\/t at fixed V^/t. For weak-to- 
intermediate coupling ADWs are present in the bulk of the trap. The agreement between the 
BALD A and the DMRG results is excellent for \U\/t < 1. With increasing \U\/t the BALD A 
scheme deteriorates 2l|, leading to an overestimation of the amplitude of the ADWs (see 
panel C). According to DMRG, the bulk ADWs disappear in the extreme strong-coupling 
limit (see panel D). For \U\/t ^> 1 a flat region of doubly-occupied sites emerges at the 
trap center, resembling that described above for the case of weak interactions and strong 
confinement (see the bottom panel of Fig. [TJ. 

The disappearance of the ADWs at strong, coupling can be explained by mapping the 
Hamiltonian (Q) onto a spin-1/2 XXZ model p, Q|, 



with J x = J y = —J z = — t 2 /\U\ and B{ = Vi(i — L/2) 2 . The total site occupation operator 
•hi is related to of by hi = 1 + of. Particle-number conservation requires working in a 
sector with fixed total magnetization (J^of) = N — L = M. In the limit \U\/t — > 00, 
J a is negligibly small and thus finding the ground state of 7ioo is equivalent to solving the 
problem of orienting a collection of spins in a nonuniform magnetic field in order to minimize 
the Zeeman energy in producing a magnetization M. Thus, for \U\/t — > 00 one expects a 
classical state with (of) = 1 ((hi) = 2) in N/2 sites at the trap center where Bi is small, 
and (af) = — 1 ((hi) = 0) in the remaining L — N/2 sites. 

Observability of the ADWs — ADWs can be detected by a measurement of the elastic 
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FIG. 1: (color online) Top panel: DMRG results for the site occupation m as a function of site 
position i for a system with N = 30 fermions in L = 180 lattice sites, and in the presence of 
a harmonic potential with Vijt = 10 -5 . For this value of Vz/t BALDA overestimates the ADW 
amplitude. Bottom panel: DMRG results (crosses) for N = 30, L = 100, and U/t = —2 are 
compared with BALDA data (filled symbols). V^/t is increased from 10 -3 to 10 . The thin solid 
lines are just a guide for the eye. 
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FIG. 2: Site occupation function of % for N = 30, L = 100, and V 2 /t = 4 X 10~ 3 . Panels 

A and B: DMRG results (crosses) are compared with BALDA data (filled circles). Panels C and 
D: DMRG results for the Hamiltonian Q (filled circles) are compared with DMRG results for the 
strong-coupling Hamiltonian © (crosses). The thin solid lines are just a guide for the eye. 

contribution to the light-scattering diffraction pattern, i.e. the Fraunhofer structure factor 

^) = ^lE e " lfc% /' (4) 

j 

through the appearance of a peak at k = k ADW . 

In Fig. 01 we show the most favorable situation where there is a wide region in the trap 
with oscillations at k ADW = ir, inducing a well denned peak in S e \(k). The height of the peak 
is non-monotonic as a function of \U \/t, as a consequence of the aforementioned crossover 
between the U — > and U — ► — oo limits, and is best observed when this ratio is of order 
two. However, in a case such as that shown in the top panel of Fig. Q where the density 
oscillations extend into regions far away from the trap center, the Fraunhofer structure factor 
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FIG. 3: DMRG data for S eX (k) as a function of k in the range 3vr/4 < k < 5vr/4, for N = 30, 
L = 100, and V 2 /t = 4 x 10 -3 . The long vertical arrow indicates the value of fc ADW . In the inset 
we show 5 e i(7r) function of \U\/t. 



peaks at a slightly lower wavenumber, /c pea k — 0.96/c ADW . 

The main results of this work, i.e. negative pair binding energies and delocalized dimers 
forming ADWs in weak traps, are not a result of finite-size effects. We have carefully 
studied the scaling of these results in the thermodynamic limit following the procedure 



proposed in Ref. 



2J, i.e. N -> oo and V 2 /t -> with N^/V 2 /t = const. We find that 
in this limit (i) Ep approaches a finite negative value and (ii) the amplitude of the ADWs 
(calculated after subtracting the smooth Thomas-Fermi site-occupation profile) approaches 
a finite value. For example, for U/t = -2 and N 2 V 2 /t = 3.6 we find E P (N > I) ft = 
-0.184 - 0.171 exp (-JV/28.61) and A ADW (N > 1) = 0.032 + 0.114exp (-iV/39.93). 

In summary, we have shown how the Luther-Emery phase emerges in an ultracold Fermi 
gas subject to parabolic trapping inside a ID optical lattice. We have demonstrated that 
the interplay between attractive interactions and harmonic confinement leads to coexistence 
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of spin pairing and atomic-density waves. The existence of a finite pairing gap A = —Ep/2 
can be tested via radio- frequency spectroscopy j^J. Atomic-density waves can be detected 
from the elastic contribution to the light-scattering diffraction pattern. At strong coupling 
atomic-density waves change into a state in which spin-singlet dimers form an extended 
region of doubly occupied sites at the center of the trap. 

This work was partially supported by FAPESP and CNPq. M.P. gratefully acknowledges 
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After this work was completed and submitted, we became aware of a related work of F. 

n 

Karim Pour et al. 24 . 



Electronic address: m.polini@sns.it 



[1] J.I. Cirac and P. Zoller, Science 301, 176 (2003); I. Bloch, Physics World 17, 25 (2004); 
D. Jacksch and P. Zoller, Ann. Phys. 315, 52 (2005); W. Hofstetter, Philos. Mag. 86, 1891 
(2006); M. Lewenstein et al, cond-mat /0606771 , 



[2] M. Greiner et al, Nature 415, 39 (2002). 

[3] T. Stoferle et al, Phys. Rev. Lett. 96, 030401 (2006); J.K. Chin et al. , |cond-mat/0607004| 
[4] B. Paredes et al, Nature 429, 277 (2004); T. Kinoshita, T. Wenger, and D.S. Weiss, Science 

305, 1125 (2004). 
[5] H. Moritz et al, Phys. Rev. Lett. 94, 210401 (2005). 

[6] T. Giamarchi, Quantum Physics in One Dimension (Clarendon Press, Oxford, 2004). 
[7] A. Luther and V.J. Emery, Phys. Rev. Lett. 33, 589 (1974). 

[8] I.V. Tokatly, Phys. Rev. Lett. 93, 090405 (2004); J.N. Fuchs, A. Recati, and W. Zwerger, ibid. 

93, 090408 (2004); A. Recati, J.N. Fuchs, and W. Zwerger, Phys. Rev. A 71, 033630 (2005). 

A Luther-Emery point has also been identified in a ID model of fermions with a narrow 

Feshbach resonance that allows them to bind reversibly into bosonic molecules: see R. Citro 

and E. Orignac, Phys. Rev. Lett. 95, 130402 (2005). 
[9] A. Seidel and D.-H. Lee, Phys. Rev. Lett. 93, 046401 (2004) and Phys. Rev. B 71, 045113 

(2005). 

[10] D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998); G.C. Batrouni et al, ibid. 89, 117203 



9 



(2002). 

[11] M. Rigol et ai, Phys. Rev. Lett. 91, 130403 (2003); M. Rigol and A. Muramatsu, Phys. Rev. 
A 69, 053612 (2004); X.-J. Liu, P.D. Drummond, and H. Hu, Phys. Rev. Lett. 94, 136406 
(2005). 



[12 
[13 
[14 
[15 
[16 

[17; 

[18 
[19 
[20 
[21 



V.L. Campo, Jr. and K. Capelle, Phys. Rev. A 72, 061602(R) (2005). 
Gao Xianlong et ai, Phys. Rev. B 73, 165120 (2006). 

E. H. Lieb and F.Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 
U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 

N.A. Lima et ai, Phys. Rev. Lett. 90, 146402 (2003). 

R.J. Magyar and K. Burke, Phys. Rev. A 70, 032508 (2004); Y.E. Kim and A.L. Zubarev, 

Phys. Rev. A 70, 033612 (2004). 

M.F. Silva et ai, Phys. Rev. B 71, 125130 (2005). 

Gao Xianlong et ai, Phys. Rev. B 73, 161103(R) (2006). 

F. Marsiglio, Phys. Rev. B 55, 575 (1997). 

This is different from the case U > 0, where the BALDA profiles agree with quantum Monte 



Carlo and DMRG data even for U/t > 1 



13J, [18|]. Since the U < Bethe Ansatz is obtained 
from that at U > by means of an exact transformation, this shows that the locality assump- 
tion inherent in the LDA is less well satisfied for attractive than for repulsive interactions. In 
fact, it appears that the BALDA performance for a fixed value of \U\/t depends on the value 
of V 2 /t 

[22] K. Damle et ai, Europhys. Lett. 36, 7 (1996). 

[23] C. Chin et ai, Nature 305, 1128 (2004); J. Kinnunen, M. Rodriguez, and P. Torma, ibid. 

305, 1131 (2005). 
[24] F. Karim Pour et ai, |cond-mat/0608491 



10 



